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All the known cases of the A/ = 4 bifurcation in the SD bands are anal- 
ysed in the framework of the microscopic theory suggested by the author 
\ [Phys. Rev. C 55 (1997) 1275]. It is shown that the intruder orbitals are of 

^ \ paramount importance to the energy staggering. 
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This talk consists of three parts: a review of the experimental situation; a short 
description of the theory; and its comparison with available experimental data. 

1. One of the recent exciting discoveries in the superdeformed (SD) rotational 
bands was the regular staggering pattern in the transition energies of the yrast SD 
band in ^'^^Gd The states of the A/ = 2 rotational band are split into two 
^ , A/ = 4 subbands separated by about 100 eV. Since then such bifurcation (or the 

\ AI = 2 staggering) has been found in the SD bands of A=130 and A=190 

^ ' mass regions. There is also some tentative evidence of the staggering in the SD 

Q ■ bands of ^^^Er |Q and ^^"^Tl In the beginning, the phenomenon showed some 



symptoms of "Pathological Science." The effect was of a magnitude close to the 
limit of detectability and had a very low statistical significance. Fantastic theories 
contrary to experience have been suggested. There were claims of great accuracy. 

Recently new measurements have been made with higher statistics using the 
Gammasphere multidetector array P, ^. The last work confirmed the original ^^^Gd 
^ . long regular staggering pattern. Besides the AI = 4 bifurcation has been observed 

^ ' in two new SD bands of ^^^Eu and ^^^Gd, which are identical to the yrast SD band of 

^■^^Gd. Both of these bands exhibit the same long and regular staggering pattern as 
^ ' the reference band. The length, regularity and striking similarity of the staggering 

■ pattern in these three bands refute the band crossing scenario [§. 

It seems natural to explain the phenomenon by the static 144 deformation of a 



nuclear shape as suggested by Hamamoto, Mottelson p and Macchiavelli et al. |jTO|. 



However, an attempt to go farther on this way and understand the staggering from 



first principles had failed In another approach, suggested first in Ref. [0], the 
AI = 2 staggering was considered to be a result of the C4 symmetry perturbation. 
The phenomenological treatment accomplished in Ref. [|l^ confirmed the dynamic 



origin of the C4 perturbation. The extension of this idea on a microscopic level was 
made in Refs. [1^, [1^ . These works are based on the mean field approach with the 
perturbation method, developed in the theory of nonadiabatic corrections to nuclear 
rotational spectra [15-17]. 

2. The starting point of the theory is the second quantized Routhian 

H' = '^Eikafak -t^Yj XaQa-Qa - ^^i-Ji - ^^2^2, (1) 

i,k '^A=2,4 
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with the quadrupole and hexadecapole residual interactions of strengths X2 and X4 
respectively. The dot stands for the scalar product Qa-Qa = J2iiQ\fj,Q\fi of the 
multipole moment operators Qx^ = J2i,k'lik{^f^)cit^k- The first term in Eq. (P is 
the single-particle energy. The last two are cranking terms involving rotation around 
an arbitrary axis perpendicular to the symmetry axis 3. Js = J2i,kjik{s)afak and 
Us, s=l,2, are the projections of the angular momentum and velocity of a nucleus. 
For simplicity we ignore the pairing interaction. 

Since all the numerical calculations done so far show the absence of any noticeable 
triaxial superdeformation, it is reasonable to confine ourselves to the axial case. 
Accordingly we perform the canonical transformation to the new single particle 
functions |z/), which satisfy the self-consistent Schrodinger equation 

XxaxoQxo - - \v) = ty\v), (2) 

A=2,4 

where j± = ji ± zj2, cij± = cui ± iu^-, and q;ao = Yl,u (Iwi/^^^ny is the axial collec- 
tive coordinate depending on the nucleon occupation number Uy. In the Hartree 
approximation, the energy in the rotating frame is given by 

r = ^e^n^ + ^^Xa^ao - ^ H XAaA^aA,-/.- (3) 

V ^ A ^ A,^j^O 

The last term of this equation contains the nonaxial collective coordinates 

«Am = Tr(gAMp), /i = ±2,±4, (4) 

which are nonzero because the density matrix p in the representation of Eq. (0) is 
not axial due to the cranking term. Since the nucleons in all orbitals except high-j 
ones are only slightly disturbed by rotation, let us decompose the nonaxial collective 
coordinates into two parts 

"Am = Qam + "Am^ (5) 

where is the multipole moment of the nucleons in the high-j shells. Being 
intruder the states of a such unique-parity shell mix weakly with other single-particle 
states. This means that a filled intruder shell does not make a contribution to 
because the operator qx^^ has a vanishing track. Therefore, is the multipole 
moment of the nucleons in the unfilled unique-parity shells near the Fermi surface. 

The second part, a\^, is the multipole moment of nucleons in the low-j orbitals. 
It is possible to use the perturbation theory to obtain the cj-dependence of this 
macroscopic quantity. Since the contribution of the intruder orbitals into the total 
nuclear multipole moment is small, Qa^j ~ ^A/^M, we have: 

«Am ^ Tr(gAMP^')) + Tr(gAy')) + .-, (6) 

where the nth order correction to the nuclear density matrix due to the 

cranking term of Eq. (j^). In performing corresponding calculations, we will use 
the representation of nonrotating axially-deformed nuclei. It is important to note 
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that the perturbation series for the axial collective coordinates, a\o, involves the 
w-independent quantity a^^o (^^^^ multipole moment of nonrotating nucleus): 



a 



(0) 
AO 



Tr(gAoP^°)). (7) 



The decomposition and Eq. (^) allow us to obtain the expansion of the 
energy (|^) in the powers of the rotational frequencies and then in the powers of 
the angular momentum components /± = Ji ± il2 (see for details Refs. fl^). In 
order to find a condition for a given band to stagger, we consider the limit of a 
pure collective rotation. Then the Hamiltonian describing nonaiabatic effects in an 
isolated rotational band has, to the quartic terms in the angular momentum, the 
form: 

H,s = AI^ + Bl^ + d{ll + J!) + c(4 + It), (8) 

The expansion (||) contains axial and nonaxial terms. The axial ones include the 
moment of inertia ^3 = 1/253) and the next inertial parameter B (see Ref. 
The nonaxial terms with the parameters 



A=2,4 



4CJ2 



involve the quadrupole and hexadecapole coupling of the intruder states with the 
nonaxial distortion (due to the rotation) of the nuclear self-consistent field. The 
latter is described by the w-independent values 

(n) (n) —(n+u)/2 —{n—ii)/2rr\ i (n\\ /'in\ 

av ="A,-M = ^+ "^'^ Tr(gA^P^ 0- (10) 

Let us note that the parameters c and d are the functions of J/Q" because of the 
tu-dependence of the values Q\^. The quadratic and quartic nonaxial terms become 
comparable for I 1^, where 

The effective rotational Hamiltonian (H) is obtained by a SU{2) mapping from 
the original fermion space to the rotor one. This method allows us to separate the 
rotational and the single-particle motion at the cost of nonadiabatic terms. For small 
spin /, when the values are close to zero, Eq. (^) is reduced to the standard 
rotational Hamiltonian of axially- symmetric nuclei, the power series of the operator 
I^. This Hamiltonian does not lead to staggering. If a nucleus has a stable nonaxial 
deformation, the collective coordinates ^22, 0:42, and 044 have cj- independent parts 
by analogy with a^Q* (see Eq.(|^)). The nonaxial terms of such rotational Hamilto- 
nian are A times greater than the ones of axial nuclei. A similar Hamiltonian has 



been considered in Ref. ||T^, where it has been shown that it generates the irregular 



A/=2 staggering. The suitable microscopic model has been studied in Ref. |T9 
For the 144 deformation with respect to the symmetry axis, ^44 is the only coor- 
dinate that involves an tu-independent part. Accordingly, the Hamiltonian has the 
C^y symmetric form of Ref. jH] . In obtaining the Hamiltonian (Rf) we neglected the 
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centrifugal effects due to collective vibrations since corresponding terms are pro- 
portional, not to Qxn, but to the coherent sums involving only a small fraction of 
intruder orbitals ||16| . 

The Hamiltonian (|^) is not fourfold invariant. The term with the operator /^+/^ 
violates the C^y symmetry. However, it is invariant under the D2 point symmetry 
group and under the transformation d —d with the simultaneous rotation by 
90° around the symmetry-axis. The later invariance allows us to consider only 
the negative values of rf, in which case the yrast band corresponds to the rotation 
around the 1-axis if / < Ic- The energy Eo{I) of the yrast states is determined in 
the quasiclassical approximation (see Ref. [Q) by the Bohr quantization condition: 

/■27r 

S^{E) = / h{(j),E)d(j) = 27Tl, (12) 





where the 0-conjugate variable Ji is the angular momentum about the axis of quan- 
tization, i.e., the 1-axis. There are two sets of the equivalent classical trajectories 
surrounding the 1-axis on the rotation energy surface. This means that the energy 
levels of the band occur in degenerate pairs. The angular momentum tunneling 
splits the degeneracy leading to the doublets. Only a single full symmetric state 
of the doublet is appropriate for the yrast SD band. Its energy is given by 

E{I) = Eoil) - 2 I T I cos{ReS2{Eo)}, (13) 

T^(^)^e-<-....(B)^/«,.E).,. (14) 

is the tunneling amplitude for the angular momentum precessing on one trajectory 
to leak over on another. The tunneling path 7 connects the points on the equivalent 
trajectories where I2 = and forms a part of the great circle, which crosses the 
saddle point at the 3-axis. Since the angle is around the 2-axis the representation 
of the Hamiltonian (H) should be used with the 2 the axis of quantization. 

Transforming the Hamiltonian to new axes we can estimate 5*2 (-E) for the classical 
minimum -Emin = {A + 2d)P + {B + 2c)I^. The action 5*2 calculated in this way 
acquires a real part if c > and Iq < I < 1^ where 



Jo = (2 - V2f/H, 



(2- v^). 



d 



nl/2 



4c 



(15) 



As spin / increases in this region, the second term of Eq. ([13]) changes a sign 
every time the ReS2 value increases by vr. Thus, we have an irregular staggering. 
For the negative value of the parameter c, the staggering is absent in the lower 
levels of the I-multiplets (the yrast band), but it exists in upper levels because the 
transformation c — — c results in the inversion of multiplet levels. It should be noted 
that the staggering is independent of the sign of c for the 6*4^ invariant rotational 
Hamiltonian. 

3. The above analysis shows that the staggering phenomenon is connected with 
the wave function oscillations in the classically forbidden regions of the rotational 
motion phase space. This effect is unusual for the second order Schrodinger equation 
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and is the characteristic of a fourth order wave equation. That is why the parameter 
c plays an important role. It involves the microscopic and macroscopic factors, Qu 
and a4^. In order to calculate the former we use the eigenf unctions of Eq. 
which, after performing the rotation of the intrinsic coordinate system by the angle 



1] around the 3-axis (cosr/ = Wi/u, ui = yuf + 0O2), has the following form for an 
intruder j-shell: 

{^nij + f^23l - i^ii^jt - (6j^ - 5)^3] - ujji} \jkr) = ejkr \ jkr), (16) 

^^^^ 2% 23 A 

^^^7(JTI)^^^' ^^^(.-1).(. + !)(. + 2)^^^' 

where Snij is the energy of the j-orbitals in a spherical potential, P2 and are the 
quadrupole and hexadecapole deformation parameters. The states of the j-shell are 
labeled hj k = 1,2, + 1/2 and the signature r = exp(— iyra). The nonaxial 
hexadecapole moment of intruder orbitals is calculated with the help of the formula: 

Q44 =510^'" I ^44 I jkr)njkr, (17) 

k,r 

where the occupation numbers of high-j intruder orbitals near the Fermi 

surface: ^13/2 for protons and ji5/2 for neutrons in the 150 and 190 mass regions. 
Figure 1 shows the dependence of this value on the number A/" of occupied orbitals 
in the j=15/2 shell and the rotational frequency u. 

Similar dependencies are observed for the j=ll/2 and 13/2 shells. The function 
Fj{M) = Q44/(n/|r'^|?7,/) is generic for the intruder configurations of the SD bands 
providing Af is small or large. This can be understood if we consider the limit of the 
angular momentum j alignment, for which the expectation value of ^44 is propor- 
tional to the Clebsch-Gordan coefficient (jm, 40|jm), with m being the projection 
of j on the rotation axis 1. The coefficient is positive for m = j and negative for 
m > j — 1, besides Z]m=i/2 0'^' "^^b"^) ~ ^■ 

The second factor, a44'', has a fixed sign. In the axially-deformed oscillator 
potential with the frequencies Uz along the symmetry axis and u± in the plane 
perpendicular to it, we have in the limit cu^ — cj^ ^ cu^ + cu^: 



(4) / 35 9h^ fuj_i + ujz'^^ 



ati = \h^ rAM2 4 2 — - 12Sx. + 3Sx^ - Saa) , (18) 

V 327r b4:M^ujj_uj^ \uj± — uJz' 

where M is nucleon mass and the quantities Sjfc are the sums of the bilinear com- 
binations of the oscillator quantum numbers rij^, ra^, A over occupied orbitals. The 
quasiclassical approximation gives us: = 3/4^2^, = 1/2^2^, Saa = 1/4^22 
for SD bands and = 3X^2, S_|_2 = Saa = for ND bands. Thus the value 
(|l^) is negative and the sign of c depends on that of Q44. 

The parameter d is calculated analogously. Both parameters c and d allow us to 
estimate the minimal spin /q, for which the AJ = 2 staggering may be present. We 
have found for the A =150 mass region the following values: Jq ~400 for the SD 
nuclei without pairing, Jq ~40 for the SD nuclei with pairing, and Jq ~10 for the 
ND nuclei with pairing. 
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We consider now a more realistic situation when the intrinsic configuration of 
a SD band has the dominating single-particle component 7r6™z/7" with m and n 
being the number of protons and neutrons occupying intruder states. By assuming 



according to Baranger and Kumar |£T[ that the mean square radius and the defor- 
mation are the same for protons and neutrons (this means that the ratio of protons 
to neutrons is the same inside a nucleus), we find for a nucleus with Z protons and 
N neutrons 



X4 



(19) 

If N is equal to Z, the last formula transforms to the parameter c of Eq. (|^). For 
the intruder configuration n6"^i>7^, the sign of c is determined by the value 

^(m, n) = 15Fi3/2(m) + 19Fi5/2(n), (20) 

which depends also on the deformation and angular frequency. The inequality JF > 
is the necessary condition for the A/ = 2 staggering to occur in the band. 

Let us apply the criterion obtained above to the explanation of experimental 
data. The simplest case is the yrast SD band of ^^^Tb, ^^^Tb(l), having the configu- 
ration 7r6^z/7^. The neutron and proton moments (^44 are negative and the staggering 
is definitely absent The values Q44(vr) and Qm{^) have different signs in the case 
of the band ^^^Dy(l) with the configuration ttG^'z/?^. This configuration leads to the 
large negative JF- value (see Table I). That is why the staggering is not observed in 
this band 

The dominant configuration 7r6^z/7^ was assigned to a number of SD bands in 
i48,i49Qj-^ and ^"^^Eu. As it follows from Table I, the neutron contribution to the 
value JF(2, 1) compensates for that of protons. In this case we have to take into 
account the contributions from other unfilled shells. Besides, the approximation of 
an isolated j-shell (|TBp is inapplicable. Nevertheless we can explain the experimen- 
tal results of Ref. qualitatively by considering ^^^Gd(l) as the reference band 
with the positive value of c. Then the identical band ^^'^Er(l) with the configura- 
tion ^''^Gd(l)®(7r[301]l/2, a = —1/2)^^ exhibits the staggering because the moment 
QaaItt) changes little after removing a proton from the Si/2 shell. The same holds 
for the band ^^^Gd{A) with the configuration i^9Gd(l)®(z/[411]l/2, a = -1/2)"^ (we 
use the labeling of the ^^^Gd bands according to Ref. ||22| and their interpretation 



according to Ref. |^). On the other hand, no statistically significant staggering has 
been observed in the bands 1 and 5 of ^^^Gd, which differ from the ^^^Gd(l) one by a 
neutron hole in the states [651] 1/2 and [642] 5/2 with the signature a = +1/2. Being 
favored these states make a noticeable contribution to the moment (^44(1^). Thus, 
the removal of these neutrons may change the sign of c. As to the band ^■*^Gd(6) 
having the configuration ^^9Gd(l)(g)(z/[651]l/2)-2 ® (z/[770]l/2, a = 1/2) at low ro- 
tational frequency and ^^^Gd(l)®(z^[770]l/2, a = —1/2)^^ at high one, the absence 
of the staggering is explained in the approximation of the dominant configurations 
7r62z/72 and nG^ur^. 
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TABLE I. Values T{2,1) for the dominant configurations 7r62z/7^( ^^^'^''^Gd) and 
7r6^z^7^ (^^^Dy(l)). The deformation parameters are taken from Ref. p^ . 



uj,MeV 








0.4 


-0.1087 


-0.1158 


-0.4190 


0.6 


0.0038 


-0.0179 


-0.2571 


0.8 


0.1454 


0.1132 


-0.3290 



There are some difficulties with the apphcation of obtained criterion to ^^^Hg 
and the cerium isotopes. Bands 1, 2 and 3 of ^^^Hg have the same intruder configu- 
ration, 7r6^i/7^. Accordingly, JF(4, 4) has a large positive value for all frequencies and 
representative deformations. It is not clear why the staggering is observed in bands 
2 and 3 and is not observed in the yrast band. The case of the cerium isotopes 
is more complicated. Having the intruder configuration 7r5'^i/6^ and J^(4, 1) > 
the yrast band of ^^^Ce shows a staggering effect. The yrast band of ^'^^Ce and 
bands 1,2 of ^'^^Ce have the configuration 7r5^z/6^ and J^{4:, 2) < for all frequencies. 
Nevertheless the staggering has been observed in Ref. 0. This observation is in a 
contradiction with the necessary condition obtained above. 

In summary, the cause of the AJ = 2 staggering has been found. The stag- 
gering is induced by the hexadecapole couphng of neutrons and protons occupying 
intruder orbitals with the nonaxial distortion of the nuclear mean field by rotation. 
By confining oneself to the pure collective rotation of nucleus we have found the 
necessary condition for the staggering in a band. All the evidences of staggered and 
nonstaggered bands in the ^4=150 mass region do not contradict with this condi- 
tion. For the configuration 7r6^z/7^, the staggering is sensitive to the occupation of 
nonintruder orbitals. 

The author is grateful to G. de France for giving the experimental results on the 
SD bands in ^"^^Gd prior to publication. This work is supported in part by Russian 
Fund for Fundamental Research through Grant No. 96-02-16115. 
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FIGURE CAPTION 

Fig. 1. Expectation value of the multipole moment Q44 in the different single-particle 
configurations of the cranked single j=15/2 shell for different rotational frequencies u!= 
0.4 MeV (circles), 0.6 MeV (squares ), 0.8 MeV (triangles), and the representative defor- 
mations /?2=0.622, /?4=0.041. 
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